Drastic facilitation of the onset of global chaos in a periodically driven Hamiltonian 
system due to an extremum in the dependence of eigenfrequency on energy 
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The Chirikov resonance-overlap criterion predicts the onset of global chaos if nonlinear resonances 
overlap in energy, which is conventionally assumed to require a non-small magnitude of perturbation. 
We show that, for a time-periodic perturbation, the onset of global chaos may occur at unusually 
small magnitudes of perturbation if the unperturbed system possesses more than one separatrix. 
The relevant scenario is the combination of the overlap in the phase space between resonances of 
the same order and their overlap in energy with chaotic layers associated with separatrices of the 
unperturbed system. One of the most important manifestations of this effect is a drastic increase 
of the energy range involved into the unbounded chaotic transport in spatially periodic system 
driven by a rather weak time-periodic force, provided the driving frequency approaches the extremal 
eigenfrequency or its harmonics. We develop the asymptotic theory and verify it in simulations. 
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A weak perturbation of a Hamiltonian system causes 
the onset of chaotic layers in the phase space around sep- 
aratrices of the unperturbed system and/or separatrices 
surrounding nonlinear resonances generated by the per- 
turbation itself [1-3]: the system may be transported 
along the layer random-like. Chaotic transport plays an 
important role in many physical phenomena [3]. If the 
perturbation is weak enough, then the layers are thin 
and such chaos is called local [1-3] . As the magnitude of 
the perturbation increases, the widths of the layers grow 
and the layers corresponding to adjacent separatrices (ei- 
ther to those of the unperturbed system, or to separatri- 
ces surrounding different nonlinear resonances) reconnect 
at some, typically non-small, critical value of the mag- 
nitude, which conventionally marks the onset of global 
chaos [1-3] i.e. chaos in a large range of the phase space 
which increases with a further increase of the magnitude 
of perturbation. Such scenario often correlates with the 
overlap in energy between neighbouring resonances cal- 
culated independently in the resonant approximations of 
the corresponding orders: this constitutes the celebrated 
empirical Chirikov resonance-overlap criterion [1-3]. 

But the Chirikov criterion may fail in time-periodically 
perturbed zero-dispersion (ZD) systems [4-6] (cf. also 
studies of related maps [7,8]) i.e. systems in which the 
frequency of eigcnoscillation possesses an extremum (typ- 
ically, a local maximum or minimum) as a function of its 
energy. In such systems, there are typically two reso- 
nances of one and the same order [9], and their overlap 
in energy does not result in the onset of global chaos 
[5-8] . Even their overlap in phase space results typically, 
instead of the onset of global chaos, only in reconnection 
of the thin chaotic layers associated with the resonances 
while, as the amplitude of the time-periodic perturbation 
grows further, the layers separate again (with a different 



topology [5-8]) despite the growth of the width of the 
overall relevant range of energy. 

The major idea of the present work is to limit the en- 
ergy range relevant to the overlap of resonances of the 
same order by chaotic layers relating either to resonances 
of a different order [10] or to separatrices of the un- 
perturbed system (these scenarios immediately suggest 
the corresponding substitute for the Chirikov criterion 
in a ZD system in the range of extremum). The onset 
of global chaos occurs respectively in a broader energy 
range or at unusually small magnitudes of the perturba- 
tion. The former case will be considered elsewhere while 
we present below details of the latter case. It relates to 
many important systems and the above effect is mani- 
fested in a very strong manner: all systems with more 
than one separatrices possess the ZD property [5,6] and 
global chaos may involve the whole energy range between 
the separatrices. 

We use as an example a potential system with a peri- 
odic potential possessing two different-height barriers in 
a period (Fig. 1(a)): 
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$ = const < 1. 



The model (1) may describe e.g. a 2D electron gas 
in a magnetic field spatially periodic in one of dimen- 
sions of the gas [11]. The interest to this system arose 
due to technological advances of the last decade allow- 
ing to manufacture high-quality magnetic superlattices 
[12,13] leading to a variety of interesting behaviours of 
charge carriers in semiconductors (see [11-15] and refer- 
ences therein). The same model describes a pendulum 
on a hinge spinning about its vertical axis [16]. 
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The separatrices of the Hamiltonian system (1) in the 
p — q plane and the dependence of the frequency lo of its 
oscillation, often called eigenfrequency, on energy E = 
Ho(p,q) are shown in Figs. 1(b) and 1(c) respectively. 
As seen from the latter, co(E) is close to the extreme 
eigenfrequency Lo m = Lo(E m ) in the major part of the 
interval [eP~\e^] while sharply decreasing to zero as 

E approaches E*jp or E^\ Such features of u>(E) are 
typical for systems with two or more separatrices. 
Let the time-periodic perturbation be added: 



q = dH/dp, p — —dH/dq, 
H(p, q) = H (p, q) - hqcos(u) f t). 



(2) 



Let us first consider the conventional evolution of chaos 
in the system (2)-(l) as h grows while ojf remains fixed at 
an arbitrarily chosen value beyond the immediate vicin- 
ity of oj m and its harmonics. It is illustrated by Fig. 
2. At small h, there are two thin chaotic layers, around 
the inner and outer separatrices of the unperturbed sys- 
tem. Unbounded chaotic transport takes place only in 
the outer chaotic layer i.e. in a narrow energy range. As 
h grows, so also do the layers until, at some critical value 
h gc = h gc (ojf), they merge. This event may be consid- 
ered as the onset of global chaos: the whole range of en- 
ergies between the barriers levels then becomes involved 
in unbounded chaotic transport. Note that the states 
{/W} ee {p = 0, q = tt/2 + 2ttZ} and {0«} = {p = 0, q = 
—n/2 + 2irl} (where I — 0, ±1, ±2, ...) in the stroboscopic 
(for instants n27r/w/ with n = 0, 1, 2, ...) Poincare section 
are associated respectively with the inner and outer sad- 
dles of the unperturbed system, and necessarily belong 
to the inner and outer chaotic layers respectively. Thus, 
the necessary and sufficient condition for global chaos to 
arise in the system may be formulated as the possibil- 
ity for the system placed initially in the state {I^ '} to 
pass beyond the neighbouring "outer" states, {C^ '} and 
{O^}, i.e. for q(t ^> 2ir/ujf) to become smaller than 
— 7r/2 or larger than 37r/2. 

A diagram in the h — u)f plane, based on the above 
criterion, is shown in Fig. 3. The lower boundary of the 
shaded area represents the function hg C (uf). It has deep 
cusp-like local minima (spikes) at frequencies ujf — Loi n ' 1 
that are slightly less than the odd multiples of uj m , 



w< n > «w m (2n-l), 



1,2,... 



(3) 



,(!) 



LO m : h gc (ujP) is 



The deepest minimum occurs at lo. 
approximately 40 times smaller than in the neighbouring 
pronounced local maximum of h gc (ujf) at w/ « 1. As n 
increases, an nth minimum becomes less deep. 

The origin of the spikes becomes obvious from the 
analysis of the evolution of the Poincare section as h 
grows while ujf w lo^\ For h = 0.001, one can see in 
Fig. 4(a) four chaotic trajectories: those associated with 



the inner and outer separatrices of the unperturbed sys- 
tem [1-3] are coloured green and blue respectively, while 
the trajectories associated with the two nonlinear reso- 
nances of the 1st order [1-3] are indicated by red and 
cyan (the corresponding attractors are indicated respec- 
tively by crosses of the same colours). Examples of non- 
chaotic (often called KAM [1-3]) trajectories separating 
the chaotic ones are shown in brown. As h increases to 
h = 0.003 (Fig. 4(b)), the blue and red chaotic trajecto- 
ries mix: the resulting chaotic layer is shown in blue. As 
h increases further (see Fig. 4(c), where h = 0.00475), 
the latter layer merges with the layer associated with the 
cyan chaotic trajectory (the resulting layer is shown by 
blue) and, finally, as h increases slightly more (see Fig. 
4(d), where h = 0.0055), the latter layer merges with the 
layer associated with the green trajectory [17], thereby 
marking the onset of global chaos as defined above. 

For most of the range [E^\ E^], w(E) is close to Lo m 
[18] (cf. Fig. 1(c)), so that the overlap between the 
resonances of the same order may be described within 
the resonance approximation. Given that it plays the 
main role in the above scenario, the function h gc (ujf) may 
also be evaluated in the ranges of the spikes using the 
resonance approximation. The explicit asymptotic (for 
$ — > 0) formulae for the minima themselves, in the lowest 
order of the parameters of smallness, are as follows: 



(2n- 1)tt 
21n(4e/$) : 



n= 1,2, 



(4) 



$ 2 <1, {lOln^e/*)}- 1 < 1 



(note that io m w 7r/{2 ln(8/$)}, so that (4) docs not 
contradict (3)), 



/i ( w (n)) = JfjL_!)££_ c « 0.179, n = l,2,..., 
9 V s ' 2V21n(4e/$) 

ln(4c/$) 
ln{101n(4c/$)/(2n- 1)} 
101n(4e/$)/(2n- 1) 

where c is the root of the equation 



(5) 



« 1, 



ln (l^xli) ~ 2x(c) = c ' x{c) " v/l = 4eC ~ 2 ' (6) 

Eqs. (4)-(6) have been derived within the following 
self-contained scheme. Let u>f be close to ncu m so that, for 
most of the energy interval [E^\e^] except in the close 
vicinity of barriers levels, the nth harmonic of eigenoscil- 
lation is nearly resonant to the driving. Then a slow 
dynamics of action I and angle tp [19] may be described 
by the following auxiliary Hamiltonian (cf. [1-8]): 

H(I,tp)= / dl(nuj — uif) — nhq n cos(tp) , (7) 

Jl{E m ) 
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where / = 1(E) = (2ir) 1 § dqp is action (note also that 
dE/dl = lo [19]) and tp = nip — Uft is slow angle [1-8] 
while q n = q n (I) is the modulus of the nth Fourier har- 
monic of coordinate (see the exact definition in [5,6]). 

For a given value of h and uif, there arc typically 
two stable stationary points in a 2ir band of the I — tp 
plane, often called resonances of the nth order. Each 
resonance is surrounded by a separatrix which includes 
one (per a 2-7T band) unstable stationary point, called 
saddle. If both saddles possess the same H the separa- 
trices reconnect each other. This condition is satisfied 
along some line u)f(h) in the h — ojf plane. We need to 
find the point (u) s ,h s ) on this line such that the lowest 
along the reconnected separatrices energy Ei coincides 
with the upper energy boundary E^ of the chaotic layer 
around the lower barrier level E^ (one can prove that 
the chaotic layer around the upper barrier level is then 
necessarily overlapped by the resonance separatrices: cf. 
Fig. 4(c)). It is known [3] that, if uif ~ u> where wo 
is the eigenfrequency at the bottom of the wells, then 
S = E^ -E^ ~ hujf/ujQ. One can prove that, if $ — ► 0, 
then the values of ui s and h s are not sensitive to the exact 
value of S within this range. Thus, the values of ui s and h s 
can be obtained using only the resonant approximation. 

In the unperturbed case (h = 0), the equations of mo- 
tion (2)-(l) can be solved explicitly, and u)(E) can be 
found then explicitly too. As <j> — > 0, it reduces to 

U W= U . \ T *E = E-\, (8) 

111 ^($-A_E)($+A_E) J 

$ < 1. 

Substituting Eq. (8) for u> into Eq. (7) for H and as- 
suming that (fcn-i V2/(2n — 1) in the relevant en- 
ergy range (note that qin = due to the symmetry of 
the potential), one can derive Eqs. (4)- (6) as the solu- 
tion of the system of two algebraic equations obtained 
from the condition of the equality of H in the saddles 
and from the condition that Ei — E^ ~ hu)f /u>o, keep- 
ing in the solution only those terms which have the low- 
est order in the parameters of smallness. The solution 
confirms the original assumptions that u> w uif/(2n — 1) 
and ^n-i ~ V2/(2n — 1) in the relevant range of ener- 
gies, which proves that the theory is self-contained. Note 
also that, to logarithmic accuracy, the asymptotic for- 
mulas (4)- (6) are valid for an arbitrary symmetric po- 
tential (whose period is normalized to 27r) and they may 
be rather easily generalized for an arbitrary Hamiltonian 
system with two separatrices close to each other. 

The values of obtained from simulations for <j> = 
0.2 (see Fig. 3) nicely agree with Eq. (4). As concerns 
h gc (wi n ^), the value $ — 0.2 is too large for Eq. (5) to 
be accurate but, even so, Eq. (5) provides a very good 
estimate for h gc (uj^p) and the correct order of magnitude 



for h gc (us ) (see Fig. 3). Calculations within the same 
scheme but using exact formulas for w and q2n-i and 
solving the system of the algebraic equations numerically 
yield the values (u>P = 0.4 ± 0.004, h gc (w^ ] ) = 0.005 ± 

0.001) and (J^ = 1.243 ± 0.01, h gc (J? ] ) = 0.025 ± 
0.008), nicely agreeing with the simulations. 
It is worth mentioning two generalizations: 

1. For non- symmetric U(q), qi n ^ and, therefore, 
pronounced spikes in h gc (u)f) exist at uif w 2ncj m too. 

2. If the time-periodic driving is multiplicative, then 
the resonances become parametric (cf. [19]). Parametric 
resonances are more complicated and much less studied 
than nonlinear resonance. But, still, the main mechanism 
for the onset of global chaos remains the same, and simu- 
lations show that h gc (uf) do possess deep spikes. At the 
same time, the frequencies of the main spikes are twice as 
large as those of the corresponding spikes in the case of 
the additive driving: this is so because the characteristic 
frequencies of parametric resonance are typically doubled 
as compared with the nonlinear resonance (cf. [19]). 

Finally, we suggest two examples of physical applica- 
tions: (i) a jump-wise increase of the dc conductivity oc- 
curs due to the jump- wise increase of the range of energies 
involved in the unbounded chaotic transport of electric 
charge carriers in a magnetic super lattice (cf. [11]); (ii) 
a significant decrease of the activation energy for noise- 
induced multi-barrier escape in the presence of periodic 
driving is associated with the noise-free transport from 
the lower barrier to beyond the upper barrier (cf. [20,6]). 

It is worth noting some analogies between our work 
and the work [21] describing the so called stochastic per- 
colation in 2D Hamiltonian systems where the merging 
of internal and external chaotic zones was also relevant. 
However, both the models and the underlying mecha- 
nisms are very different (there is no zero-dispersion situ- 
ation in [21]; rather the two-dimensionality is inherently 
important); besides, in [21], there was no explicit quan- 
titative description analogous to our Eqs. (4)-(6). 

The work was supported by INTAS Grants 97-574 and 
00-00867. We are grateful to J. Howard, P. McClintock, 
J. Meiss and K. Richtcr for discussions. 
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FIG. 1. The potential U(q), the separatrices in the phase 
space and the eigenfrequency uj{E) for the unperturbed sys- 
tem (1) with $ = 0.2, in (a), (b) and (c) respectively. 
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FIG. 2. The evolution of the stroboscopic (for instants 
n2ir/ujf with n — 0, 1,2, ...) Poincare sections of the system 
(2)-(l) at <3> = 0.2, with uif = 0.3 while h grows from the 
top to the bottom: (a) 0.01, (b) 0.05, (c) 0.1. The number 
of points in each trajectory is 2000. In parts (a) and (b), 
three characteristic trajectories are shown: the inner trajec- 
tory starts from the state {/'°'} = {p = 0, q — vr/2} and is 
chaotic but bounded in the space; the outer trajectory starts 
from {O' '} = {p = 0, q = — 7r/2} and is chaotic and un- 
bounded in coordinate; the third trajectory represents an ex- 
ample of a regular (KAM) trajectory separating the above 
chaotic ones. In (c), the chaotic trajectories mix. 
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FIG. 4. The evolution of the stroboscopic Poincare section 
of the system (2)-(l) with $ = 0.2 as the amplitude of the per- 
turbation h grows while the frequency is fixed at ujf = 0.401. 
The number of points in each trajectory is 2000. The chaotic 
trajectories starting from the states {1^} and {O' '} are 
drawn in green and blue respectively. The stable stationary 
points (the lst-order nonlinear resonances) are indicated by 
the red and cyan crosses. The chaotic layers associated with 
the resonances, in those cases when they do not merge with 
those associated with the green/blue chaotic trajectories, are 
indicated in red and cyan respectively (their real width is 
much less than the width of the drawing line). Examples of 
KAM trajectories embracing the state while separating 

various chaotic trajectories are shown in brown. 
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FIG. 3. The bifurcation diagram indicating (by shading) 
the perturbation parameters ranges at which global chaos ex- 
ists. The integration time for each point of the grid is 120007T. 
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